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Abstract 

Non-equilibrium phase transitions are studied theoretically for the two-dimensional Hubbard model subject 
to bias voltages from the electrodes coupled to the system. By combining the fluctuation exchange approxi- 
mation with the Keldysh method for non-equilibrium, we have studied the properties of the non-equilibrium 
Fermi liquid phase and determined the phase diagram with transition to non-equilibrium magnetic and 
superconducting phases. 
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1. Nonequilibrium Superconductivity 
Electron Correlation 



and 



Due to the decades of research for the HTC 
cuprates, organic superconductors and other classes 
of materials, the understanding of correlated electron 
systems is becoming matured. One novel avenue to 
pursue is alternative methods for doping other than 
the conventional chemical doping. One way is to in- 
ject carriers from interfaces or electrodes which has 
been recently realized in oxides heterostructures[l|. 
If we consider the case where two electrodes are at- 
tached, then by adjusting the bias voltage across the 
electrodes, we can introduce carriers whose distribu- 
tion may be different from equilibrium. The effect of 
such non-equilibrium carriers in correlated electron 
systems has been theoretically proposed recently by 
the present authors Q • The main interest is whether 
superconductivity can occur in non-equilibrium in 
an open, two-dimensional correlated electron system 
coupled to electrodes (Fig. [1] (a)). We have adopted 
the Hubbard model as a standard model for corre- 
lated electron systems. The non-equilibrium situa- 
tion is treated with the Keldysh formalism, while the 



many-body effect with the FLEX approach^- We 
also note that non-equilibrium methods for doping 
cuprates have been experimentally realized in ref. [3( 
where photoinduced superconductivity in HTC ma- 
terials was observed. Our approach is expected to 
shed light to these systems as well. 

2. Nonequilibrium Phase Diagram and Super- 
conducting Order 

We consider a thin layer of strongly correlated 
material described by the two-dimensional Hubbard 
model which is coupled to electrodes. Here we assume 
that the electrodes are attached in the vertical direc- 
tion (Fig. Q] (a)) and double step Fermi distribution 
is realized in the correlated system (Fig. [5]). 

In non-equilibrium systems, an important effect of 
the electron-electron interaction is the smearing of 
the electron distribution function. As can be seen in 
Fig. [5J the double step Fermi distribution f® s real- 
ized by the electrode and bias V becomes smeared 
into a smooth non-equilibrium distribution function 
f c s. In this non- equilibrium Fermi liquid phase the 
bias voltage V effectively plays a similar role as the 



Preprint submitted to Physica C 



December 7, 2010 




Figure 1: (a) Two dimensional correlated system (shaded) in 
a finite bias V. (b) Obtained non-equilibrium phase diagram 
against the bias (V), temperature (T), and the doping level 
(<5) with antiferromagnetic (AF) and d-wave superconducting 
(SC) phases. 

temperature since the smearing of the distribution 
becomes stronger as one increase V. 

The non-equilibrium phase diagram obtained 
within the Keldysh + FLEX method is schemati- 
cally shown in Fig. [1] (b). Away from half- filling, 
for which FLEX is expected to be applicable, the su- 
perconducting gap function with d-wave symmetry 
dominates even in finite bias voltages. However, as 
one increases the bias, we can see that both the AF 
and SC phases become suppressed. Indeed, this can 
be seen in Fig. |3] (b). First at bias V = 0.1, which is 
near the non-equilibrium Fermi liquid to d-wave su- 
perconductor transition, the AF fluctuation is sharp 
and has multiple peak structure which reflects the 
incommensurate nesting. However, as one goes far 
away from the transition (V = 0.2), the peak be- 
comes smaller and the structure vanishes. 

The suppressed orderings can be understood in 
terms of the change in the non-equilibrium elec- 
tron distribution. The AF fluctuation is induced by 
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Figure 2: Nonequilibrium distribution function for two values 
of the bias V at half filling (S = 0). Dashed lines are the 
noninteracting distribution function f® s . 

electron-electron scattering across the Fermi surface 
which is enhanced near half-filling due to nesting ef- 
fect. However, in the presence of finite bias V, (i) the 
Fermi surface becomes split into two with Fermi ener- 
gies Ep ± V/2 where Ep is the original Fermi energy, 
and (ii) smearing takes place. The first effect reduces 
the nesting while the second weakens the scattering. 
The AF fluctuation is the glue for the transition from 
(non-equilibrium) Fermi liquid to d-wave supercon- 
ductivity and the reduction of the fluctuation results 
in the destruction of the superconducting phase. 

In conclusion, we have studied the transition be- 
tween the non-equilibrium Fermi liquid, AF magnetic 
order and d-wave superconductor which is control- 
lable by the applied bias. Finally we comment on the 
smearing effect which destroys the long range order in 
finite bias. In the Keldysh + FLEX method, the self- 
consistent treatment of the Green's function tends to 
overestimate the smearing effect. If the smearing is 
not as strong, there is a possibility that a new order 
appears which does not exist in equilibrium 
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Figure 3: (a) Superconducting gap function Re<p(k,Lj = 0) 
in non-equilibrium. The bias is V = 0.1, 0.2. Dashed lines 
represent nodes, (b) Spin susceptibility x[Qi w = 0) f° r ^ = 
0.1, 0.2. For S = 0.14, U = 4.5, and fi = -0.35. 
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